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ABSTRACT 

Motivated by coronal mass ejection studies, we construct general relativistic 
models of a magnetar magnetosphere endowed with strong magnetic fields. The 
equilibrium states of the stationary, axisymmetric magnetic fields in the magne- 
tar magnetosphere are obtained as solutions of the Grad-Shafranov equation in a 
Schwarzschild spacetime. To understand the magnetic energy buildup in the mag- 
netar magnetosphere, a generalized magnetic virial theorem in the Schwarzschild 
metric is newly derived. We carefully address the question whether the mag- 
netar magnetospheric magnetic field can build up sufficient magnetic energy to 
account for the work required to open up the magnetic field during magnetar 
giant flares. We point out the importance of the Aly-Sturrock constraint, which 
has been widely studied in solar corona mass ejections, as a reference state in un- 
derstanding magnetar energy storage processes. We examine how the magnetic 
field can possess enough energy to overcome the Aly-Sturrock energy constraint 
and open up. In particular, general relativistic (GR) effects on the Aly-Sturrock 
energy constraint in the Schwarzschild spacetime are carefully investigated. It is 
found that, for magnetar outbursts, the Aly-Sturrock constraint is more strin- 
gent, i.e., the Aly-Sturrock energy threshold is enhanced due to the GR effects. 
In addition, neutron stars with greater mass have a higher Aly-Sturrock energy 
threshold and are more difficult to erupt. This indicates that magnetars are prob- 
ably not neutron stars with extreme mass. For a typical neutron star with mass 
of 1 — 2M©, we further explore the effects of cross- field current effects, caused 
by the mass loading, on the possibility of stored magnetic field energy exceeding 
the Aly-Sturrock threshold. 
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1. INTRODUCTION 

After the discovery of soft gamma repeaters and anomalous X-ray pulsars (Mazets et al. 
1979; Mereghetti & Stella 1995), magnetar models of these sources are proposed to explain 
the relevant phenomena (Duncan & Thompson 1992; Thompson, Lyutikov & Kulkarni 2002). 
Magnetars are believed to be neutron stars with strong magnetic field, ~ 10 14 — 10 15 G 
(Duncan & Thompson 1992). The magnetar outbursts, such as giant flares, occur with huge 
release of magnetic energy ~ 10 44 — 10 46 ergs. The energy for magnetar outbursts is widely 
accepted to be supplied by the star's magnetic field. However the physical process by which 
the energy is stored and released is one of the great puzzles in high energy astrophysics. Two 
possibilities exist for the location where the magnetic energy is stored prior to an eruption: 
in the magnetar crust or in the magnetosphere. For the former possibility, a giant flare may 
be caused by a sudden untwisting of the magnetar interior magnetic field (Thompson & 
Duncan 2001). Subsequently, a sudden and brittle fracture of the crust leads to the giant 
flare. In this crust scenario, the energy stored in the external twist is limited by the tensile 
strength of the crust. Alternatively, based on the short timescale of the giant flare rise time, 
~ 0.25ms (Palmer et al. 2005), the second possibility — the magnetospheric storage model, 
was proposed by Lyutikov (2006). The energy released during an eruption is stored slowly 
(on a longer timescale than the timescale of giant flare) in the magnetar magnetosphere prior 
to the eruption. An abrupt reconfiguration and dissipation of the magnetic field due to a loss 
of confinement (Flyer et al. 2004) or a dynamical instability (Lyutikov 2003; Komissarov et 
al. 2007) produces the giant flare. This mechanism has the feature that the energy stored 
in the external twist may not be limited by the tensile strength of the crust, but instead by 
the total external magnetic field energy. 

The magnetospheric storage model of magnetar giant flare shares similar magnetic en- 
ergy buildup process to solar eruptions, such as coronal mass ejections (CMEs). In this 
model, the energy released during an eruption is stored in the magnetospheric magnetic field 
before the eruption. Large-scale eruptive CMEs often give rise to the opening up of magnetic 
field lines that were originally closed. The processes of magnetic fields opening up have been 
extensively investigated in the CME studies (Barnes & Sturrock 1972; Aly 1984; Mikic & 
Linker 1994). It is physically reasonable to assume that the preemption closed state must 
possess more magnetic energy than the posteruption open state. As will be discussed in 
detail below, requiring the magnetic field to open imposes an extreme energy constraint on 
theories for CMEs. This energy requirement on solar CMEs has been under extensive theo- 
retical studies in the past decades (Aly 1984; Sturrock 1991; Wolfson & Dlamini 1997; Zhang 
& Low 2005). The energy storage processes take place quasi-statically on a long timescale. 
When the magnetic field reaches a threshold, due to the instability or loss of confinement, 
the field erupts suddenly on a much shorter dynamical timescale. Analogous processes of 
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magnetic field opening up are believed to occur in magnetar giant flares (Woods et al. 2001; 
Thompson et al. 2002; Beloborodov 2009). All these features of the storage model are in 
good agreement with the observations of magnetar giant flares (Lyutikov 2006). 

The similarity between solar eruptions and magnetar giant flares (Lyutikov 2003) mo- 
tivates this study on the energy buildup process in the magnetar magnetosphere. We note 
that there are important differences between solar eruptions and magnetar outbursts. For 
situations in the magnetar magnetosphere (Beloborodov & Thompson 2007), the location 
where the magnetic energy buildup occurs is quite near the neutron star surface (~ 1 — 2i? NS ). 
General relativistic (GR) effects near the neutron star surface are important (Ciolfi et al. 
2009). General relativistic effects are currently, however, not taken into account in rele- 
vant energy storage processes. In this work we will investigate these processes with GR 
spacetime curvature effects considered. More specifically, we will ignore effects of magnetar 
rotation since they are slow rotators and describe the background geometry of the magnetar 
magnetosphere with the Schwarzschild metric. 

The virial theorem is a helpful tool for us to understand the energy properties in the 
magnetar magnetosphere. The flat spacetime magnetic virial theorem (Chandrasekhar 1961) 
has been extensively exploited in astrophysical researches (Aly 1984; Zhang & Low 2005). 
Attempts to get the GR virial theorem have been made by Chandrasekhar (1967), but he 
just considered a hydrostatic system, with the effects of magnetic fields completely ignored. 
In this study we establish the magnetic virial theorem in the Schwarzschild metric, which 
helps us to better understand GR effects on the physical behaviors of the magnetic energy 
buildup. 

For the magnetospheric storage model, an important question for the giant flare ener- 
getics is: Can the magnetospheric magnetic field store enough energy before an eruption? 
For the magnetic energy alone to power a magnetar giant flare, the energy must be sufficient 
to open up the magnetic field. However, for the nearly force-free magnetic field exterior to 
a sphere, a well-known result by Aly (1984, 1991) and Sturrock (1991) suggests that the 
energy of a fully open field is the upper limit on the energies of all the force-free fields in 
simple geometries^]. Thus the transition from a closed field configuration to an open one 
(which is actually required for a realistic eruption) is not energetically favored. Due to 
this Aly-Sturrock constraint, the initial magnetic field before eruption must have energy in 
excess of the threshold set by the Aly-Sturrock energy constraint. This Aly-Sturrock con- 
straint is widely discussed in the solar CMEs studies. But its implications for magnetars 



1 Here simple geometries mean that the two ends of all field lines are anchored onto the neutron star 
surface. 
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are only briefly mentioned in Lyutkov (2006). Furthermore, GR effects on this important 
Aly-Sturrock constraint have not been considered in prior works. One purpose of this work 
is to clarify how GR effects influence the Aly-Sturrock constraint. 

The Aly-Sturrock constraint constitutes a bottleneck for the storage model of magnetar 
giant flares. There are a number of ways, however, to avoid the Aly-Sturrock constraint. A 
deviation from a perfectly force-free initial state might make a difference. In this scenario, it 
is expected that the cross-field electric currents are viable source of energy for the eruption. 
Detailed calculations about solar CMEs by Low & Smith (1993) suggested that a non-force- 
free magnetic field with cross field currents due to the mass loading of plasma can store 
more energy than the Aly-Sturrock field. Such mass loading effects are further discussed 
by Wolfson & Dlamini (1997) and Zhang & Low (2004). The mass loading of plasma in 
a non-force-free magnetic field acts like a rigid wall to confine the magnetic field, in other 
words, it would act as a lid that allows the magnetic energy to increase above the limit, 
and when the lid is suddenly removed, the field springs outward (Fan & Low 2003). By 
analogy, it is possible that in the magnetar magnetosphere, the mass loading plays the same 
role to compress the magnetic field. Consequently, the magnetic field can store magnetic 
energy above the Aly-Sturrock constraint. But no theoretical calculations were performed to 
corroborate this idea. In this work we will provide such a demonstration. Another possibility 
for the magnetic energy to exceed the Aly-Sturrock constraint is the formation of detached 
field lines from the magnetar surface (magnetic bubble or magnetic flux rope, e.g., Low & 
Smith 1993; Flyer et al. 2004), which will be further discussed in Yu et al. (2011, in prep). 

This paper is organized as follows: in §2 we introduce the generalized magnetic virial 
theorem in the Schwarzschild metric. In §3 we will discuss how the the Aly-Sturrock field 
energy is affected by general relativistic effects. We will explore the cross-field effects caused 
by the mass loading on the magnetic energy storage in §4. Conclusions and discussions are 
given in §5. 

2. Generalized Virial Theorem in Schwarzschild Spacetime 

The virial theorem is of vital importance for understanding the magnetic energy storage 
in the magnetar magnetosphere. In the flat spacetime, it was proposed by Chandrasekhar 
(1961) and has been used widely in solar physics researches (e.g., Low & Smith 1993). We 
focus in this paper on the physical behavior near the magnetar surface, GR effects should 
be incorporated. Because observed magnetars have a very slow rotation rate, we ignore 
the rotation effects and adopt the Schwarzschild metric as the background spacetime. In 
this section we establish the virial theorem in the Schwarzschild metric including effects of 
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magnetic fields. We consider a steady state magnetosphere around magnetars. The metric 
g^y of Schwarzschild geometry reads (Misner, Thorne & Wheeler 1973) 

ds 2 = g^dx^dx u = -a 2 dt 2 + a~ 2 dr 2 + r 2 d6 2 + r 2 sin 2 6d<p 2 . (1) 

The factor of a is defined as 

a(r) = sjl- 2 -^ , (2) 

where r g = GAi QS /c 2 is the gravitational radius, G is the gravitational constant, Ai ns is the 
mass of the neutron star, and c is the speed of light. 

A plasma containing only a perfect fluid and an electromagnetic field, is described by 
the energy-momentum tensor (Anile 1989) 

= TZ d + T£m = (p + P + b 2 ) u»u» + (p + |) «T - , (3) 

where p is the isotropic pressure, p = po + ppry is the energy density (including that due 
to the rest mass p ) and b 2 = b^. A polytropic equation of state is adopted and we take 
7 = 4/3 throughout this paper. Here the Einstein summation rule is assumed and Greek 
letters take on the values t, r, 9, and 0. The magnetic field 4-vector is 

h n = * F ^ Uv ? ( 4 ) 

where *F fJiU is the Maxwell tensor and u v is the four velocity of the comoving observer 
(Anton et al. 2006). The plasma is assumed to be in magnetostatic equilibrium, thus the 
four velocity is = ((— gtt) 1 ^ 2 , 0, 0, 0). Under such circumstances, the condition 



V„T^ = , (5) 



reduces to 



*r-V^ + (p+p o +6 V ^p^ + (v=g w) + T*yr = o , (6) 

y dx v 2 y dx v y/=gdx v vv ; Xa ' u 

where g is the determinant of the metric g^ u and the explicit expressions of the connection 
coefficients r^ CT (Weinberg 1972) are given in Appendix A. 

A re-arrangement of various terms of the above equation using Gauss theorem leads to 
the following generalized virial theorem in the Schwarzschild spacetime (Details are given in 
Appendix A), 



E + (37 — A)U = J a 2 (^-+pj ( r ' rfS ) 
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a\B ■ r)(B-dS) + / (1 f ] ( B 2 + B 2 + ^-Jip ) dV . (7) 



Jv 2 V 7 7-1 
In this equation, r is the position vector. Here dS is a surface area element directed outwards 
and dV is a volume element, both measured by a locally inertial observer. The factor of a 
is given in equation (j2J). Note that the total energy E is the sum of the magnetic, internal 
and gravitational potential energy, namely 

E = M + U + W, (8) 

where 

M = J T dV ' (9) 

U^J-^dV, (10) 

w = -j PoG ^ dV, (11) 

are the magnetic, internal and gravitational potential energy, respectively. Here we have 
absorbed a 47r factor into the definition of the magnetic fields throughout this paper. In the 
above equations, the magnetic field B in the "ordinary" orthogonal basis (defined in Section 
4) is used. The relation between B and the magnetic field 4-vector 6 M is given explicitly 
in Appendix A. Note that B r is the radial component of B and B 2 = B 2 + B 2 + B 2 . 
Throughout this paper, we mainly work with the magnetic field B. This choice is made 
mainly for the convenience of comparison between the results in the curved spacetime and 
the flat spacetime. 

The last integral on the right hand side in equation (J7J) appears owing to general rel- 
ativists effects. This term disappears when taking the flat spacetime limit, i.e., a 2 — > 1. 
Note also that this equation becomes the usual virial theorem in the flat spacetime as a 2 — > 
1 (Chandrasekhar 1961). In particulai@, for the magnetically dominated force-free field, we 
arrive at 

M= f ^ (r-dS) - [ a 2 (B • r) (B-dS) + / (1 ~ ^ (B 2 + B 2 ) dV . (12) 

Assuming that the magnetic field vanishes sufficiently rapidly at large distances, we find that 
the energy of the force-free fields in the exterior r > r of the neutron star is 

M = tttq [ a 2 (B 2 - Bf[ sin OdO + / ^ (2B 2 . + B%) dV , (13) 



2 Although we will be treating non- force- free magnetosphere in which cross-field effects (caused by mass 
loading) are important, the discussion here is restricted to magnetically dominated force-free fields. The 
relevance will become clear as we proceed. 
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where tq is the radius of the neutron star. 

We note that, in a flat spactime, the second term on the right hand side of the above 
equation disappears and the total magnetic energy of a force-free magnetic field in the 
exterior region r > tq of a sphere is uniquely determined by the field values at the the 
boundary r = tq. However this is no longer the case for the curved spacetime, since additional 
terms proportional to r g appear on the right hand side of this equation. Close observation of 
equation ( 1T5j) shows that, when GR effects are ignored, no force-free field that is completely 
detached from the solar surface (i.e., B r = at r = r in the exterior region r > r ) can exist 
(Low 2001). However, such completely detached field configurations in the general relativistic 
magnetar magnetosphere, due to the spacetime curvature, may be in the equilibrium statd^l 
This suggests that, besides the normal flux at the magnetar surface, the GR spacetime 
curvature provides additional self-confining effects. As a result, it needs more work to be 
done to open the magnetic field in the curved spacetime than in the flat spacetime. It is 
conceivable that when the magnetar mass increases, this effect becomes more evident (see 
Figure [2]). Such GR effects have important implications for the magnetic energy storage 
process in the magnetar magnetosphere. In the next section, we will quantitatively calculate 
their influences on the Aly-Sturrock constraint. 



3. Aly-Sturrock Constraint For Magnetic Field Energy 

To discuss the magnetic energy in the magnetar magnetosphere, it is beneficial to in- 
troduce the potential field B pot in the Schwarzschild metric which satisfies (Uzdensky 2004) 

V x (aB) = , (14) 

and the boundary condition 

r = r ,B r = F{6) . (15) 

In this paper we mainly discuss the dipole field and its relevant open state. The explicit 
expression of the dipole field can be found in Appendix B. In this case the above boundary 
becomes B r — C cos 6, where C is a constant. Note that the potential field now involves the 
spactime curvature term a in equation (j!4p . This is quite different from the flat spacetime 
definition (Komissarov 2004). Note that, as a— >1, this potential field definition reduces to 
its flat spacetime form. The associated magnetic energy of the potential field is designated 
as M pot . For the force- free field in the magnetosphere, there exists one interesting energy 



3 See the magnetic field configuration in Figure 8b of Low (2001), which can not maintain equilibrium in 
flat spacetime. But such configurations can be self-confined by the spacetime curvature effects. 
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reference state, the Aly-Sturrock state (Aly 1984,1991; Sturrock 1991). Imagine all force-free 
magnetic fields complying with the boundary condition f JT5|) . with one end of each line of 
force anchored to the star's surface and the other out to infinity. Among all these fields, the 
one with the lowest energy is potential everywhere except for a current sheet at the equator 
(Aly 1984,1991; Sturrock 1991). This lowest energy state is the Aly-Sturrock state. Call this 
magnetic field configuration B open . The total energy of this state is designated as M open . 
The well-known Aly-Sturrock conjecture claims that for any fully closed force-free fielcEl with 
the boundary condition (IT5|) . its total energy M FF satisfies the following relation, 

M pot < M FF < M open . (16) 

This first half of this inequality means a current-free potential field is the lowest energy 
state. And the second half suggests that the opening up process of an initial closed force- 
free magnetic field requires considerable amount of work to be done on the magnetic field. Of 
particular interest is whether the pre-eruption magnetic energy M can exceed the threshold 
set by the Aly-Sturrock field. This is crucial for the magnetically driven outbursts. 

Some numerical experiments have recently demonstrated the validity of this conjecture 
(Antiochos, DeVore & Klimchuk 1999; Hu 2004). Due to the importance of the Aly-Sturrock 
constraint for the magnetic eruption, it is worthwhile to reconsider this problem when GR 
effects are important. Note that this Aly-Sturrock state is unique (Aly 1984; Sturrock 1991) 
and can be constructed by the following technique. Modify the boundary condition (1151) to 

r = r ,B r = \F(9)\ . (17) 

After getting the field with this boundary condition and reversing the directions of those 
lines at the boundary r = r$ where B r < of this field, we could get the Aly-Sturrock state 
(see also Low & Smith 1993). We have calculated the fully open field B open and the relevant 
energy M open numerically. The details to obtain the Aly-Sturrock field and the magnetic 
energy M open are discussed in Appendix C. In Figure [TJ, an illustrative example of the fully 
open Aly-Sturrock field is shown. The current sheet at the equator is shown by a thick solid 
line. 



4 Strictly speaking, this condition is not fulfilled since the field lines open at the light cylinder. Fortunately, 
magnetars are slow rotators, so the light cylinder is quite far away from the neutron star surface and this 
effects can be negligible. For this reason, we focus in this paper on the non-rotating neutron stars. 
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Fig. 1. — The fully open Aly-Sturrock field with an initial dipole field boundary condition. 
The thick solid line at the equator denotes the current sheet in the field. 
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3.1. Dependence of M open on Neutron Star Masses 

To investigate the spacetime curvature effects on the Aly-Sturrock constraint, we calcu- 
late the Aly-Sturrock threshold M open for different magnetar masses. Throughout this paper 
we take the neutron star radius r = 1, so for a neutron star mass of 1 — 3 M & , r g ranges 
from 0.15 — 0.45 (For simplicity, we keep the neutron star radius fixed at 10 km, though this 
is not the case in reality). In Figure [21 we show the variation of M open (in units of M pot ) with 
the neutron star mass. This figure shows that the more massive the magnetar, the higher 
the threshold is. For instance, for the dipole field with r g = 0.15 (1 M Q ), the energy of the 
fully open Aly-Sturrock field is M open = 1.80M pot ; when r g = 0.21 (1.4 M ), the energy 
becomes M open = 1.88M pot . Consequently, it is more difficult for more massive neutron stars 
to surpass the Aly-Sturrock energy threshold. From this figure, we also note that as r g — > 
the Aly-Sturrock threshold approaches the flat spacetime limit M open = 1.662M pot . This is 
consistent with our physical expectation. 

This increase with mass of the Aly-Sturrock energy threshold stems entirely from the 
spacetime curvature self-confining effects mentioned in Section 2 and this behavior is quite 
different from the solar eruption in flat spacetime, in which the Aly-Sturrock field energy (~ 
1.662 M po t) is independent of the star mass. It should be emphasized that in the magnetar 
outbursts, the Aly-Sturrock energy constraint is more stringent than for the solar CME-type 
eruptions. From Figure 2, we can infer that magnetars are probably not neutron stars with 
extreme mass ~ 3M , as the Aly-Sturrock threshold could hardly be reached. 

For typical neutron star masses (~ 1 — 2M , r g ~ 0.15 — 0.3), it is necessary to seek 
initial magnetic fields which possess magnetic energy in excess of the threshold set by the 
Aly-Sturrock energy M open . One possibility is due to the mass loading effects. The estimated 
ejected mass loading is about 10 22 g (Lyutikov 2006). This mass loading can be balanced by 
pressure forces in the vertical direction. The pressure gradient in the horizontal direction, 
however, requires magnetic forces associated with cross-field currents, i.e., J x B ^ 0, to 
maintain the equilibrium state. The deviations from a strictly force-free magnetic fields, i.e., 
the cross-field contribution, are worth further investigations (Low & Smith 1993, Wolfson & 
Dlamini 1997). Physically speaking, the mass loading would act as a lid over the magnetic 
field. The field can be compressed globally by a sufficient amount of plasma. As a result, the 
energy of the compressed magnetic field increases as the total load increases, and eventually 
the magnetic energy exceeds M open . In other words, the cross field current densities provide 
additional sources of magnetic free energy which may be enough to enable the magnetic field 
to clear the threshold M ODen . 
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Fig. 2. — The variation of the M open (in units of M pot ) with r g (equivalently, the magnetar 
mass). Note r g ranges from 0.15 to 0.45, which corresponds to a mass range of 1 — 3M . 
Note that in flat spacetime the fully open field energy M open = 1.662M pot for the dipole field 
boundary condition, which is denoted by dot-dashed line in this figure. 
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4. Axisymmetric Magnetostatic Magnetosphere with Cross-Field Currents 

In this section we explore the cross-field effects on the magnetic energy storage proper- 
ties in the magnetar magnetosphere. Similar investigations in solar CMEs have been carried 
out by Zhang & Low (2004). Specifically, we will focus on the question whether the mag- 
netic energy in the magnetosphere can exceed the Aly-Sturrock threshold. In what follows, 
we consider that the magnetar magnetosphere evolves quasi-statically on sufficiently slow 
timescale that we can treat the magnetosphere as being essentially in magnetostatic equi- 
librium. A steady state axisymmetric, purely poloidal magnetic field in the Schwarzschild 
metric can be written as 

B = B pol = W x V0 , (18) 

where \I/(r, 9) is the poloidal magnetic stream function. The "ordinary" orthogonal basis is 
used, where e^ = (Weinberg 1972, no summation rule over \i is used in this equation), 

namely, 

1 1 

e ? = ad r , e § = ~d e , e; = — — . (19) 
r v r sin V 

The poloidal magnetic field components are (Uzdensky 2004) 

B = —n -1S7T> ""IT- • ( 20 ) 



r sin 6 \r 36 ' dr J 

To account for the the cross-field current effects induced by the mass loading, we must go 
beyond the force- free approximations (Yu 2011) and turn to the full magnetohydrodynamic 
(MHD) equation (J6j). This equation decomposes into the following two equations 

d ( 2 <9*\ sin# d ( 1 <9*\ 2 2n dp , s 

Wr\ a ^) + —w\^re^e) +r sm ^ = ' (21) 

rr dp GAlns . , 

9 fr + ( p + ti—^ = ' (22) 

for balance across and along the magnetic field (Low & Smith 1993). A simple solution to 
equation (1221) reads 



1 P(*) 

Po ~GM n „ r™ 



m I I 2m - 2 - 7 > Tg 



7 — 1 J r 



(24) 



where P(*) is a free function of the magnetic stream function * and m is a constant. 



To keep the problem mathematically tractable, we take the free function P to be linear 
in *. Subsequently, equation (123]) and (124"]) become 

P=^. (25) 



13 



Po 



1 A(^ + ^ ) 



GM, 



ra+l- 2m + 2 + 



7 



(26) 



7 — 1 J r 

where \l/ an d A are constants. Substitute equation fl25|) into equation f l2~Tj) . we obtain the 
following linear Grad-Shafranov equation 



8 ( 2 8^\ sind 8 
dr \ dr J r 2 86 \ sin 6 86 



1 <9#\ sin 2 fl 

+ A— — r = 



-m— 1 



The general solution to the above equation can be written as 

^ = f m (r) sin 2 6 + ^ pot , 
where ^ po t is an arbitrary potential stream function satisfying (Ghosh 2000) 

sin^ 8 ( 1 d^ po t 



8_ / 2 ^pot 
dr \ dr 



+ 



r 2 86 V sin 6 86 







(27) 



(2? 



(29) 



This equation can be readily solved by the variable separation method (see Appendix B). 
The function f m (r) satisfies the following second order ordinary differential equation (ODE) 



\ tv* J /y 2 } ' ~ 



A 



.m-l 







(30) 



where prime denotes derivatives with respect to r. The particular solutions can be readily 
obtained analytically. For m = 3,4, 5, 6, 7, and 8, the radial function f m are given explicitly 
in Appendix D. The simple linear solutions given by the above equations can not be expected 
to describe the magnetar magnetosphere in realistic details. However, the solution given by 
equation (T2"g|) with ^/ pot = can be used to obtain a physical estimate of how much energy 
can be stored in the magnetosphere prior to eruptions. 

The magnetic energies for different values of m and r g are listed in Table 1. We find that, 
for r g = 0.15 and 0.21 , the field configurations are able to sustain magnetic energy higher 
than the Aly-Sturrock threshold as m > 8. Simple estimation shows that the total magnetic 
energy of a magnetar with magnetic field ~ 10 14 — 10 15 G is approximately 10 46 — 10 48 ergs. 
Given the actual giant flare energy release ~ 10 44 — 10 46 ergs, we know that a few percent 
of the magnetic energy in excess of the Aly-Sturrock threshold is needed to release during 
a giant flare. This energy requirement can be fulfilled as m reaches a critical value. For 
instance, we find that for r g = 0.15 the magnetic energy M with m = 8 is approximately 15 
percent above the Aly-Sturrock threshold M open , which is enough to drive magnetar giant 
flares. In Figure 3, we show the m = 8 solution with r g = 0.15 



fs(r) 
fs(r ) 



sin 2 6 , po 



1 A^ 

GM n .~r™ 



m + 1 



2m + 2 + 



7 



7 
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where the stream flux is normalized to unity at r = r$ and 6 = 71/ 2 and we have set the 
constant \l/o — in equation fl26|) . The left panel in this figure shows the magnetic field lines 
and the right one shows the contour of the density departure^] from an arbitrary, spherically 
symmetric distribution. In this particular state, the magnetic energy M is 2.12M pot , which 
is greater than the Aly-Sturock state M open = 1.80M pot . 

The solution with m = 3 and \l/ pot = is a purely radial magnetic field, 

B r = A ^ 2 , Bq = . 

This solution has been extensively discussed in the Blandford-Znajek process (e.g. Blandford 
& Znajek 1977) related to relativistic astrophysical jets. But in our discussion this state itself 
is of no particular interest as we are more concerned with the initial closed state. To introduce 
the closed field structures, we add a dipole field to the m = 3 purely radial magnetic field, 
i.e., 

* = A^-sin 2 ^±^ dipole , (31) 
Wo) 

where the stream function is also normalized. The magnetic fields with r g = 0.15 are shown 
in Figure 4. The left panel in this figure corresponds to the "+" sign, which approximately 
models the effects of the neutron star wind (Bucciantini et al. 2006). Such configurations 
are also discussed by Low & Tsinganos (1986) and applied to model the effects of solar wind. 
Note that when A increases to 5.0, the magnetic energy in the left panel is about 1.83M pot , 
exceeding the corresponding Aly-Sturock energy M open by about 2%, which suggests this 
state may support a giant flare. If A is even increased, more magnetic energy can be obtained. 
The right panel takes the "— " sign in the above equation. Though the right panel shows a 
state that is physically unacceptable, it is worth pointing out that the energy of the state 
with detached field lines (~ 2.86M pot ) is much higher than the energy in the left panel. This 
also suggests that, when there are detached fields in the magnetosphere, the stored magnetic 
energy can be much larger than those configurations whose field lines are all anchored to the 
magnetar surface. This possibility to bypass the Aly-Sturrock constraint has been discussed 
by Flyer et al. (2004) for solar CMEs and will be further discussed for magnetar giant flares 
(Yu et al. in prep). 



5 Note that an arbitrary, spherically symmetric density distribution corresponds to the term that is pro- 
portional to the constant "Jo in equation (|26j) . which is ignored in this hgure. 
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X X 



Fig. 3. — The m = 8 solution with r g = 0.15. The left panel shows the magnetic field. The 
right panel shows the density departure from an arbitrary, spherically symmetric distribution, 
i.e., we set the constant \l/ = in equation (j26|) . The energy of this state is 2.12M pot , which 
exceeds the Aly-Sturrock threshold 1.80M pot . 
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X X 



Fig. 4.— The magnetic fields in equation ( 1311) with A = 5.0 and r g = 0.15. The left 
panel corresponds to the "+" sign and the right the "— " sign. The energies of the two 
state are 1.83 and 2.86 M pot , respecitvely. Both of them exceed the Aly-Sturrock threshold 
M opcn = 1.80M pot . Though the state in the right panel is physically unacceptable, it is 
instructive that configurations with detached field lines can build up more energy than the 
simple connected field lines. 
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5. Conclusions and Discussions 

We construct general relativistic models of non-rotating neutron stars endowed with 
strong magnetic fields. The equilibrium states of axisymmetric force-free magnetic fields 
in magnetar magnetospheres are found as solutions of the Grad-Shafranov equations in the 
Schwarzschild geometry. A newly derived general relativistic magnetic virial theorem in 
presented in this work. Based on this magnetic virial theorem, we carefully examine the GR 
effects on the well known Aly-Sturrock energy threshold. We found that this energy threshold 
increases with the magnetar mass. As a result, it is more difficult for massive magnetars 
to erupt. By this observation, we conclude that magnetars are probably not neutron stars 
with extreme mass. The non-force-free magnetic field induced by the mass loading is further 
investigated as a possibility to bypass the Aly-Sturrock constraint for typical magnetar mass 
around ~ 1.4M Q . 

We mainly discuss dipolar surface boundary conditions in this paper. This is the case 
for magnetar's large scale fields. But observations show a striking feature that the emergence 
of a strong four-peaked pattern in the light curve of the 1998 August 27 event from SGR 
1900+14, which was shown in data from the Ulysses and Beppo-SAX gamma-ray detectors 
(Feroci et al. 2001). These remarkable data may imply that the geometry of the magnetic 
field was quite complicated in regions close to the star where GR effects are important. As 
a result, complex boundary conditions should be important for the outburst of magnetars. 
Effects of different boundary conditions on the energy buildup in magnetars are worth further 
investigations (Antiochos et al. 1999). 

For simplicity, we have neglected the relativistic wind from the neutron star surface. 
Actually, the wind from the neutron star (e.g. Bucciantini et al. 2006) may cause part of 
the magnetic field lines to be in the open states before eruption. Similar effects have been 
explored in solar CMEs (Low & Smith 1993; Wolfson 1993). It is interesting to investigate 
the effects of neutron star wind on the magnetic energy storage properties. 

Helicity has been discussed extensively in solar physics (Zhang & Low 2005). CMEs 
are believed to be the unavoidable products of the coronal evolution as a result of magnetic 
helicity accumulation (Zhang et al. 2006). But helicity in the GR regime is not a well- 
explored issue. Finding a self-consistent definition of helicity in the curved spactime and 
investigating the relevant helicity properties are interesting topics for further explorations. 

The field topology change from a closed state to an open state must be accompanied by 
the magnetic reconnection. After a certain threshold is reached, the dynamical instability 
sets in. The gradual quasi-static evolution of the magnetar's magnetosphere will be replaced 
by the dynamical evolution of the field. This naturally explains the problem as to how 
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a very slow buildup of the external shear (over an interval of ~100 yr) could lead to the 
sudden release of external magnetic energy on a much shorter timescale (Lyutikov 2006). 
The magnetic energy dissipation in the strongly magnetized plasma is caused by the tearing 
mode instability (Lyutikov 2003, Komissarov et al. 2007). Relativistic tearing instability 
induced reconnections in the nonlinear regime need further studies to better understand the 
magnetar outburst behaviors. 

Our theoretical models can not address the nonlinear dissipation processes that occur 
during giant flares. However, current GRMHD simulations provide a unique opportunity 
to study the dynamical outburst physics. The models constructed in this work are likely 
to be useful as initial states in GRMHD numerical simulations to explore the dynamics of 
magnetic eruptions (Gammie et al. 2003, Yu 2011). 

We thank the anonymous referee for important comments and suggestions that improve 
this paper greatly. The research is supported by the Natural Science Foundation of China 
(Grant 10873033, 10703012, 10778702 and 10973034), the Western Light Young Scholar 
Program and the 973 Program (Grant 2009CB824800). The computation is performed at 
HPC Center, Kunming Institute of Botany, CAS, China. 
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Table 1: Values of M/M pot for different values of m. 



m 


r g = 0.15 


r g = 0.21 


r g = 0.3 


3 


2.19 


2.30 


2.55 


4 


1.00 


1.00 


1.00 


5 


1.15 


1.13 


1.10 


6 


1.44 


1.38 


1.30 


7 


1.77 


1.67 


1.52 


8 


2.12 


1.98 


1.76 


9 


2.48 


2.29 


2.00 


10 


2.84 


2.61 


2.25 



The gravitational radius r g is taken as 0.15, 0.21 and 0.3, which correspond to magnetar 
mass of 1.0, 1.4 and 2.0 M Q . The Aly-Sturrock energy thresholds, shown in Figure 2 for the 
three values of the magnetar mass, are M open = 1.80M pot , M opcn = 1.88M pot and M open = 
2.06M pot , respectively. According to this table, we note that, when m > 8 for r g = 0.15, 0.21 
and m > 10 for r g = 0.3, the magnetic energy in the magnetosphere could be higher than 
the Aly-Sturrock threshold. 
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A. Derivation of Virial Theorem in Schwarzschild Metric 



The four equations expressing conservation of energy momentum are 

V„T^ = , 



(Al) 



where the Einstein summation rule is assumed and Greek letters take on the values t, r, 6, 
and <p. The four velocity for a plasma in magnetostaic equilibrium is 

«* = (-9a)- 1/2 , u r = u 9 = u* = 0. (A2) 

Given the energy-momentum tensor in equation (151) , the covariant derivative can be expanded 
as follows, 

W = <r^L (p + ~) + r S> + p + & V« CT 



1 d{yf=g b"b u ) 



-fj 



dx h 



(A3) 



The radial component of the above equation becomes (note that the connection coefficients 

tt 



pM _ 1 „p,v dgtt \ 
1 " _ 2^ flx"-' 



., <9_ / /r 
\ 

1 9 



(p + p + b 2 ) 



dg u 
dr 







-g dr 



(V^g VV) 



+ 



-9 d0 



(v 7 ^ b r b e ) + T r rr b r b r + T r 9e b% e + T^b^j = 



(A4) 



where 



-(r-2r s ), r^ = -(r-2r g )sin 2 ^. 



r(r — 2r 9 ) 

The "ordinary" component of the magnetic field B in the orthogonal basis (Weinberg 1972) 
is related to the magnetic field 4-vector b^ by 

B r = ^ b r = b r ,B e = yfgTe b e = b e , = Jg£ b* = ^ b^ , (A5) 



and 



b u W = B 2 = B 2 + Bl + Bl . 



(A6) 



Multiplying the equation (1A4|) by r and expressing the magnetic field 4-vector 6 M by the 
"ordinary" magnetic fields B in equation (1A4[) . we may arrive at 

d 

I n -1 I - — I n 4 n .... />'- 1 - 

Tn 



a 2 r 



Or 



(p+^-)+^(p + P + B 2 )-^l {rWBl) - ^4 (sin ^ 



sin 



+ ^Sjf + o a (^ + S^ + ^) = 0. 



(A7) 



Performing the volume integral with the usage of Gauss's theorem, the above equation can 
be re- arranged to give the generalized virial theorem, equation (j7J) in the main text. 
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B. Dipole Field Boundary Conditions and Separable Solutions for Potential 

Fields 



To get the dipole field boundary conditions, we need to obtain the current-free potential 
field. To be self-contained, we describe the separable solutions of the homogenous Grad- 
Shfranov equation, which are also the building blocks for the Aly-Sturrock fully open field. 
The homogenous GS equation reads 



d_ 

dr 



1 



2r g ^ <9$ 
dr 



+ 



sin 8 d 



1 



r 2 86 V sin 6 06 







Separable solutions of the above equation are of the form 

*(r,0) = R(r)Q(9) . 
Substitute the above equation into equation (IB 1 j) . we obtain 



d 



1 dQ 



d6 V sin 9 dQ 



-A- 



e 

sin 6 



d 
dr 



1 _2r £ \dR 



x R 



(Bl) 



(B2) 



(B3) 



(B4) 



r j dr 

where A is the separation constant. The lowest order of solution are the special case with 
A = 0, which can be obtained by setting A = in the above two equations. The solutions 
are then 

&(6) = acosfl + b , (B5) 



R(r) = c , 

where a,b, and c are constants. This solution is the Schwarzschild monopole. 



(B6) 



The order of the solution is denoted by the ordinal number m (m — 1 corresponds to 
the dipole field), related to the constant A by A = m(m + 1). Equations flB3j) and (IB4j) 
become (Ghosh 2000): 



'I 



d 2 e 

1 - V 2 )-r^ + m ( m + 1)6 = 
dfi 2 

9 d 2 R dR 
z 2 )^- - 2— + m m + 1)R -- 

dz 2 dz 



(B7) 
(B8) 



where \i = cos 9 and z = r jr g — 1 . 
The solution of equation ( 1B7h is 



00*) = (i - /i 2 



dPm(fl) 

dji 



(B9) 



-22 - 



where P m {^) is the Legendre polynomial. The solutions of equation (IB8j) are 

?£3(*) 



i2(r) 



(BIO) 



where V^\{z) and Q^|(z) are Jacobi polynomial and Jacobi functions of the second kind, 
respectively. For r 3> r g , the Jacobi polynomial and Jacobi function's asymptotic behaviors 
are (Szego 1939) 

r^l(z) ~ r™- 1 , q£3(*) ~ r— a . (Bll) 

The superscripts in the Jacobi polynomial and Jacobi function will be suppressed hereafter, 
as the values remain the same throughout this study. The explicit expressions for the Jacobi 
polynomials and Jacobi functions can be found in Gradshteyn & Ryzhik (1980). 



Of particular interest is the dipole configuration (m = 1) determined by 



:i-^ 2 



r 2 Q (z) 



r . ( r 

— In 

2 V r - 2rv 



sin 2 fl 



(B12) 



This solution can be used as boundary conditions. 



C. Determination of the Aly-Sturrock Field 

To appreciate the Aly-Sturrock constraint on the availability of free magnetic energy, 
we need to determine the Aly-Sturrock state numerically. Following Low & Smith (1993), 
the boundary conditions of Aly-Sturrock fully opened field can be obtained by flipping the 
flux function according to the boundary condition (JTj 



* / *(n>,0) o<e<n/2 

^modify - | ^ (r0) 7r/2) _ ^ Q) ^2 < < 7T ' K ' 

Specifically, for the original dipole boundary condition, the modified boundary condition 
becomes 

, n \ r, a f sin2 ^ 0<9<tt/2 ,_ n . 

*-*('»•») = * { 2 - rf( ^2 <1 <'„ • < C2 > 

where 

and ro is the magnetar radius. The solutions to the homogeneous Grad-Shafranov equation 
are of the form 



*(r, 9) = a n (r 2 Q n _a(r)) (sin 2 O^j^) + a + a, cos9 

n=l ^ A* / 



(C3) 
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where \x = cos 9 and Q n -i{r) is the Jacobi function of the second kind. It is clear that 

a = BoAt , a x = -B A X . (C4) 
We define the following flux function as 

^j*(r,9) = \&(r, 9) - a - axcosO . 
The problem becomes to determine the coefficient a n 



(C5) 



n=l ^ " 

subject to the modified boundary condition (1C2|) 

^*(r , 0) = *(r , 0) - a - ax cos 9 

D f sin 2 0-l + cos0 O<0<vr/2 

= ^ X ( l- s in 2 + cos0 ff /2<fl< W • (C6) 

According to the orthogonality of associated Legendre polynomials -P^(/i), we have that 

1 2n + 1 



f\*{r ,B)P l Mde , (C7) 
>/ o 



r 2 Q n _ 1 (r )2n(n + l)7 

Note that \l/*(ro,6 l ) is an odd function of in the integration range. When n is an odd 
integer, the coefficients a n 's vanish. The non-zero coefficients a ra 's (n = 2m) can be written 

as 

° 2m = ' rga!llVo)2m 4 p m + + 1 l) f 2 < Sm29 - 1 + ™0)PL(™0)dO , (C8) 

where P 2m (cos0) is the associated Legendre polynomial. After some manipulations, we arrive 
at 

= B Ax 4m +1 (-l) m - 1 (2m-2)! = c 2m 
tt2m riQ 2m _x(r ) m(2m + 1) 2*»(m - l)!(m + 1)! " Q 2m -i(r ) " 1 J 

The radial and 9 components of the magnetic field, according to equation (T20]) . are 

oo 

B r = «2mQ 2m _i(r) [2m(2m + l)P 2m ] - ^ , (CIO) 



m=l 

and 



Be = -yjl - ^ f> m [2Mr)+rO] sin0^|^ , (Cll) 

m=l 
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where prime denotes derivative with respect to r. The magnetic energy of the open field, 
according to equation ffTBT) . is 



M opcn = vrr 



(l M ^B^'+YV ' 



Am(2m + 1) 



2m(2m + l)-(l-^| f 2 + ?#*# 



Am + 1 

2" 



47rr 5J B 2 A? / -dr + Anr g ^ c: 

Jr m=1 



., [2m(2m + l)] z 
2m 4^TI 



Q2m-l(r) 
Q.2m-l{r ) 



r dr 



+ 2vrr, 



2 4m(2m + 1] 



'2m 



4m + 1 



(r - 2r 9 ) 



2Q 2m -i(r)+rQ 2m _ 1 (r) 
Qam-i(»"o) 



m=l -■- . - jr 

The potential dipole field energy M pot can be calculated as follows, 

B T = 2B g r (r) cos 9 , 
B e = B g e (r) sin 9 , 



dr 



where 



1 f r 

q r (r) — - In 

y V ; 2 \r-2r, 



r r 
1 9 a 



g e {r) = M 1 - 



2r„ 



2r g (r - r g ) ^ 
r(r — 2r n ) 



r — 2r, 



The potential dipole field energy is 

M pot = \J Bl {Ag 2 r (r) cos 2 9 + g 2 (r) sin 2 9) dV . 



(C12) 

(C13) 
(C14) 

(C15) 
(C16) 

(C17) 



D. Solutions for the Ordinary Differential Equation ( 1301) 

For m = 3, 4, 5, 6, 7, and 8, the functions / m (r)'s are 

J3 2 , 



/4 = A 



2rr 9 + 2r 2 + r 2 ln(r — 2r 5 ) — r 2 lnr 



8r| 



/ 5 = A 



6r 2 r 9 + 6rr 2 + 8r^ + 3r 3 ln(r — 2r g ) — 3r 3 In r 



48rr^ 



(Dl) 
(D2) 
(D3) 
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6 r 3 r„ + 6r 2 r 2 + 8rr 3 + 12r 4 + 3r 4 ln(r - 2r„) - 3r 4 In r 

/. = A ^-15^ ■ < D4 > 

30r% + 30r 3 r 2 + 40r 2 r 3 + 60rr 4 + 96r„ 5 + 15r 5 In , s 

= 2880^7 ^ ■ < D5 > 

30r 5 r„ + 30r 4 r 2 + 40r 3 r 3 + 60r 2 r 4 + 96rr a 5 + 160r a 6 + 15r 6 In 

/8 = 7680rt^ ' (D6) 

respectively. In the calculation of the magnetic energy in the exterior of the neutron star, 
the stream functions are normalized as 

q=Ur) sin 2 d (D7) 
Jm{ro) 
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